Operator perturbation theory in the backward Heisenberg picture 
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We present a simple operator perturbation theory in the backward Heisenberg picture. Compared 
with the well-known Heisenberg picture, the revised picture is based on the backward time instead 
of the forward time. The unique feature of the uncommon picture is that the perturbation expansion 
becomes very simple and the famous Dyson expansion is not directly involved. Its relationship with 
the perturbation expansion for the density operator developed by Kubo is also discussed. 
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Introduction. Consider a perturbed quantum system with the Hamiltonian 

H(t)=H + H 1 (t), (1) 

where Ho is the time-independent part and Hi is a small time-dependent perturbation which is assumed to be switched 
on at initial time 0. We arc interested in the following question: given an observable F that does not explicitly depend 
on time and its Heisenberg picture defined as 

F H (t) = U*(t)FU(t), (2) 

where U(t) is the time evolution operator of the Hamiltonian H(t) and the superscript H stands for the Heisenberg 
picture, can we expand the dynamic observable in the order of the perturbation based on the unperturbed system 
without directly resorting to the famous Dyson expansion [l|? Apparently, if the constraint was dropped, the problem 
would become almost trivial because of 

U(t) = u {t)T+eW- 1 ti^W Hl W u '>W* r (3) 

where 7+ is the time-ordering operator, and Uq{t)~ cxp(Hot / ih~) is the time-evolution operator for the unperturbed 
system. This question was raised when we investigated the fluctuation relations in nonequilibrium processes 04^] , 
and in the work we give a definite answer. We must emphasize that we aim to reinvestigate the old operator 
perturbation issue from a very different viewpoint and are not intended to present a method as a substitution for the 
Dyson expansion. In fact, we cannot even completely exclude that the same method has been in the literature. 

The backward Heisenberg picture. We start by introducing a revised Heisenberg picture observable @ 

F B (t,t') = U(t')F H (t)U\t') (4) 

with < t' < t. Because t' is the backward time, we call the picture as the backward Heisenberg picture and 
specifically indicate it by a superscript B. Obviously, if one selects t'=0 and t, the picture reduces to the Heisenberg 
and Schrodinger pictures, respectively. The equation of motion for the new defined dynamic observable F B (t, t') with 
respect to t' is simply 

ifid t ,F B (t,t') = [H{t'),F B {t,t')]. (5) 

It is worthy emphasizing that this is a terminal condition problem with F B (t, t)=F and that the Hamiltonian operator 
therein is indeed the Schrodinger picture, which is significantly different from the standard Heisenberg equation of 
motion, 

iW t F H (t) = -[H H (t),F H (t)]. (6) 

with the initial condition F H (0) = F. The careful reader may notice that Eq. ([5]) is almost same with the equation 
of motion for the density matrix p{t) except that t' here is in place of the forward time t, which is not accidental 
because of 

d e Tv{F B (t,t')p(t')}=0. (7) 
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The expectation value (F)(t) of the operator at time t in the picture is 

Tr{Fp(t)} = Tr{F H (t)p } = Tr{F B (t, 0)po}, (8) 

where po is the initial density operator that may be a pure or mixed state. To illustrate the uncommon picture, 
we calculate the position and momentum's backward Hciscnbcrg pictures for a one-dimensional harmonic oscillator 
with a Hamiltonian H = p 2 /2m + muj 2 x 2 /2. Because the system is time-independent and one may easily find their 
equations given by 

d t ,x(t,t') = -p(t,t')/m, (9) 
d v p{t,t') = muj 2 x{t,t'). (10) 

For simplicity in notations, we did not write the superscript B out explicitly. Solving both equations with the terminal 
conditions x(t, t) = x and p(t, t) = p, we obtain 

P P 

x(t,t') = \xsvaujt cos Lot] sin wr + [x cos Lot H coswtl coswi', (11) 

mui rruo 

p(t,t) = — [uilox sinwt — pcoscot] cosujt' + [mux coscot + pcostot] sin Lot . (12) 

Letting t' = 0, the above solutions reduce to those directly calculated in the Heisenberg picture @. In fact, if the 
Hamiltonian does not depend on time, the backward Heisenberg picture ((H) could not possibly show any distinctive 
features compared with the standard Heisenberg picture because of obvious 

F B (t,t') = F H (t-t'). (13) 

Operator's perturbation. Now we turn our attention to the perturbation case ([T]). Substituting the perturbation 
expansion 

F B (t, t') = F B {t, + F B {t, t') + F 2 B (t, t') + ■ ■ ■ (14) 
into the equation of motion ([5]) and comparing its both sides order by order, we find 



(15) 
(16) 



ihd t ,F B (t,t') = [H ,F B (t,t>)] 
F B (t,t) = F, 

m>F B (t,t>) = [H 0l F B (t,t')] + [^(f),^ (t,f),] 
F B (t,t) = 0, 

i%d v F B {t,t') = [Hq, F B {t, t')] + [HiV^FfW)] n7) 
F B (t,t)=0, [U) 

Their formal solutions are 

F B (t,t') = U,{t')F B {t)ul{t')=F B {t-t'), (18) 
F B (t,t') = -(iR)- 1 ^ (Oy'*c4(^)[^i(^)»^o fl (*y)]^o(T , )dT / Cr t(/), (19) 

F B (t,t') = -(ihr^t 1 ) f U$(T J )[H 1 (T'),F*(ty)]Uo(T')dT'uZtf), (20) 



respectively. The simplicity of Eqs. (fT5ll -(fl7 |) contrasts sharply with the complexity of the perturbation equations for 
the Heisenberg picture observable 

F H (t) = F B (t) + F B (t) + F B (t) + ■ ■ • , (21) 
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each term on the right hand side of which satisfies 



ihd t F H (t) = -[H ,F H (t)} 
F o H (0)=F, 



Ff(Q) = 0, 

' ihd t F 2 H (t) = 

{ F 2 H (0)=0, 



[H ,F»(t)) m^it) + (ih)- 1 fZ[H ,(H l )g(r)]dT,Fi l (t) 



(22) 
(23) 



-[^ 0!J Ff(i)]-[(i?i)^(i) + (^)- 1 /o[^o J (ffi)f(r)]dr ) Ff( i ) 



-[(^-^VoVO^rO^ (24) 
-(^)- 2 /o(^i)?(r')^ifo/o(^i)?(^)^", Ff(t)] 



respectively, where [ , ] + denotes an anticommutator. We expand them until the second order, because the evolution of 
the n-th order F^(t) depends on all orders from to n, and the calculations of these coefficients become considerably 
tedious as the order increases. However, we must point out that the relation © has imposed conditions 



F£(t) = F*(t,0) (n = 0,1,2,.-.). 



(25) 



In the Schrodinger picture, Kubo has developed a method to evaluate the perturbation expansion for the density 
matrix 



p(t) = Po + Pi(t) + p 2 (t) + 



(26) 



on the basis of its equation of motion in 1957, which does not depend on the Dyson expansion as well |8(. One may 
see that our calculation above is highly analogous to his method. Particularly, Eq. ((SJ has indicated 



Tr{p n (t)F}=Tr{F*(t,0)p o }. 



(27) 



It is interesting to concretely calculate how these identities arise using the solutions (fl~8 |) -([2"0 |) . We use the first three 
order terms (n = 0, 1, 2) as examples. The identity with zero order is obviously trivial. For the first and second order 
terms, we have 

Tr{F 1 B (t,0) Po } 
= -iiK)- 1 /''Tr{[ J ff 1 (T'),F ( f(t-T')]po}dT' 



(th)- 1 / TT{[H 1 (r'),p }F Q H (t-T')}dr' 



Tr {(ih)- 1 f U (t - r')[#i(T'), A)]*4(* ~ T 'W F h 
Jo 



(28) 



and 



Tr{F 2 B (t,0)p } 



t r t 



W- 2 J q / ^oV) [i?i(r'),C/o(r'-r")[i?i(r"),^(r"-0]t/o T ( 

" 2 f [ T Tr{[/ V) \ Hi {t")Mt" -rO^iCr'),^^ -« T ( 
Jo Jo 
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V-r") U {T')pv}dT'dT'- 



(ih) 

Tr {(ih) 



oV'-r 



JO 



U (r") Po }dT'dT" 

U (t - t') [H^t'), U (r' - (A p ]U\r' - r")] U T Q (t - r')dr'dr" F}, 



(29) 



respectively. Two simple formulas Tr{ [A, B]C} = Tr{[B, C]A} and Tr{[A, [B, C}} D} = Tr{[B, [A, D]C}} are used in 
the derivations. We immediately find that the integral terms before the operator F in the last equalities of these 
equations are just pi (t) and p 2 (t) , respectively [1, 0| . 



Conclusion, we give a simple operator perturbation theory in the backward Hcisenbcrg picture. Its extension of 
the classical mechanical systems is straightforward. 
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